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Microscopic Analysis of Order Parameters in Nuclear Quantum Phase Transitions 
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Microscopic signatures of nuclear ground-state shape phase transitions in Nd isotopes are stud- 
ied using excitation spectra and collective wave functions obtained by diagonalization of a five- 
dimensional Hamiltonian for quadrupole vibrational and rotational degrees of freedom, with pa- 
rameters determined by constrained self-consistent relativistic mean-field calculations for triaxial 
shapes. As a function of the physical control parameter - the number of nucleons, energy gaps 
between the ground state and the excited vibrational states with zero angular momentum, isomer 
shifts, and monopole transition strengths, exhibit sharp discontinuities at neutron number N = 90, 
characteristic of a first-order quantum phase transition. 

PACS numbers: 21.60.Jz, 21.60.Ev, 21.10.Re, 21.90.+1 



Phase transitions in equilibrium shapes of atomic 
nuclei correspond to first- and second-order quantum 
phase transitions (QPT) between competing ground- 
state phases induced by variation of a non-thermal con- 
trol parameter (number of nucleons) at zero tempera- 
ture. Theoretical studies have typically been based on 
phenomcnological geometric models of nuclear shapes 
and potentials, or algebraic models of nuclear struc- 
ture [U, 0, |j|, but more recently several attempts have 
been made towards a fully microscopic description of 
shape QPT startin g fr om nucleonic degrees of freedom 
H IS ISM IS B US EH- In particular, in Refs. @, H we 
have reported a microscopic study of nuclear QPT in the 
region Z = 60, 62, 64 with N ~ 90, based on constrained 
self-consistent relativistic mean-field calculations of po- 
tential energy surfaces. While in Rcf. [7] the generator 
coordinate method (GCM) was used to perform configu- 
ration mixing of angular-momentum and particle-number 
projected relativistic wave functions restricted to axial 
symmetry, in Q collective excitation spectra and tran- 
sition probabilities have been calculated starting from a 
five-dimensional Hamiltonian for quadrupole vibrational 
and rotational degrees of freedom, with parameters deter- 
mined by constrained mean-field calculations for triaxial 
shapes, i.e. including both (3 and 7 deformations. The re- 
sults reproduce available data, and show that there is an 
abrupt change of structure at TV = 90 that can be approx- 
imately characterized by the X(5) analytic solution at the 
critical point of the first-order quantum phase transition 
between spherical and axially deformed shapes. 

A phase transition is characterized by a significant vari- 
ation of one or more order parameters as functions of 
the control parameter. Even though in systems com- 
posed of a finite number of particles, i.e. in mesoscopic 
systems, phase transitions are actually smoothed out, in 
many cases clear signatures of abrupt changes of struc- 
ture properties are observed. In their study of QPT tran- 



sitions in mesoscopic systems [12j, Iachello and Zamfir 
have shown that the main features of phase transitions, 
defined for an infinite number of particles, N — > 00, per- 
sists even for moderate N ~ 10. Their analysis has 
been followed by a number of studies of shape phase 
transitional patterns in nuclei as functions of the num- 
ber of particles, e.g. number of bosons in the frame- 
work of intcracting-boson-type models. As emphasized 
in Ref. [T3|, there are two approaches to study QPT: (i) 
the method of Landau, based on potentials, and (ii) the 
direct computation of order parameters. In the case of 
atomic nuclei, however, a quantitative analysis of QPT 
must go beyond a simple study of potential energy sur- 
faces. This is because potentials or, more specifically, de- 
formation parameters that characterize potential energy 
surfaces, are not observablcs, and can only be related to 
observables by making very specific model assumptions. 
The direct computation of observables related to order 
parameters has so far been based mostly on particular nu- 
clear structure models, e.g. the interacting boson model, 
in the framework of which such observables are defined 
as expectation values of suitably chosen operators. 

In this work we combine both approaches in a con- 
sistent microscopic framework, and present an illustra- 
tive example of calculation of observables that can be 
related to quantum order parameters as functions of nu- 
cleon number. An order parameter is a measure of the 
degree of order in a system. As a normalized quantity 
that is zero in one (symmetric) phase, and non-zero in 
the other, it characterizes the onset of order at the phase 
transition. When symmetry is broken, several variables 
can be introduced, related to order parameters, to de- 
scribe the state of the system. As in our previous studies 
@, [U , the shape transition in Nd isotopes with N ^ 90 
will be considered. The analysis starts by performing 
constrained self-consistent relativistic mean-field calcu- 
lations for triaxial shapes, i.e. including both (3 and 
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FIG. 1: (Color online) (a) Self-consistent RMF+BCS triaxial quadrupole binding energy map of 150 Nd in the ^ - 7 plane 
(0 < 7 < 60°). The contours join points on the surface with the same energy (in MeV). (b) Projection of the binding energy 
map on the 7 = 0° axis, (c) The dependence of the binding energy on the deformation parameter 7, for two values of the axial 
deformation j3 = 0.2 and 0.25. 



7 deformations. The resulting self-consistent solutions, 
i.e. single-particle wave functions, occupation probabil- 
ities, and quasiparticle energies that correspond to each 
point on the binding energy surface, are used to calcu- 
late the parameters that determine the collective Hamil- 
tonian: three mass parameters, three moments of inertia, 
and the zero-point ener gy c orrections, as functions of the 
deformations f3 and 7 [1 31 ] - The diagonalization of the 
Hamiltonian yields the excitation energies and collective 
wave functions, that are used to calculate observablcs. 

In Fig. [T] we plot the self-consistent triaxial quadrupole 
binding energy map of 150 Nd in the (3 — 7 plane (0 < 
7 < 60°). The calculations are performed by imposing 
constraints on expectation values of the quadrupole mo- 
ments (Q20) an d (Q22), the relativistic functional PC-F1 
(point-coupling Lagrangian) [Til ] is used in the particle- 
hole channel, and a density-independent J-force is the 
effective interaction in the particle-particle channel, with 
pairing correlations treated in the BCS approximation. 
As shown in 0] , the binding energy maps of Nd isotopes 
show a gradual transition from lighter spherical nuclei to- 
wards the strongly prolate deformed 152 Nd. Of particular 
interest are nuclei around 150 Nd, for which experimental 
evidence for shape phase transitional behavior has been 
reported [Hj]. 150 Nd is considered to be a good example 
of empirical realization of the X(5) model for the criti- 
cal point of first-order phase transition between spherical 
and axially deformed shapes [lij . 

The microscopic binding energy surface of 150 Nd dis- 
plays a flat prolate minimum that extends in the interval 
0.2 < (3 < 0.4 of the axial deformation parameter (Fig. 
lb), and a parabolic dependence on 7 for 7 < 30° in the 
region of the flat prolate minimum (Fig. lc). The flat 
bottom of the potential has been considered a signature 
of possible phase transition. There are, however, prob- 
lems when one considers deformations as possible order 
parameters of phase transitions (lo| . First, deformation 
parameters are not observablcs, and can only be linked 



to observables, e.g. transition rates, excitation energies, 
within the framework of a specific model. Second, in cal- 
culations that include both (3 and 7 degrees of freedom, a 
phase transition cannot be characterized by the behavior 
of just one deformation parameter, e.g. (3. Results ob- 
tained with the five-dimensional collective Hamiltonian 
show that many properties of the excitation spectra are 
affected by (3 — 7 coupling. Band-head excitation ener- 
gies, energy spacings within the bands, and transition 
strengths depend on the 7 stiffness of the potential [l7j ]. 

Here we analyze obervablcs that arc directly computed 
using collective wave functions obtained from a micro- 
scopic five-dimensional collective Hamiltonian. The im- 
portant question is how much are the discontinuities at 
a phase transitional point smoothed out in finite nuclei 
and, second, how precisely can a point of phase transi- 
tion be associated with a particular isotope, considering 
that the control parameter, i.e. nucleon number, is not 
continuous but takes only discrete integer values. Fig. [5] 
displays the differences between squares of ground-state 
charge radii (r 2 ) +(A + 2) — (r 2 ) +(A), and the isomer 

shift (r 2 ) 2 + — (r 2 ) + between the first 2 + state and the 
ground state, as functions of the neutron number. The 
former displays a peak at N = 88, 90, whereas a pro- 
nounced discontinuity is predicted for the latter between 
N = 88 and N = 92. 

Signatures of ground-state phase transitions in quan- 
tum systems characterize the evolution of both excitation 
spectra and order parameters. In Refs. 18, 0| the 
scaling properties of the energy gap between the ground 
state and the first excited vibrational states with zero 
angular momentum was studied for a system of Nb in- 
teracting bosons. At the critical point of the phase tran- 
sition the gap is strongly reduced in finite systems, and 
goes to zero as N — > 00. In the left panel of Fig.[3]we plot 
the isotopic dependence of the first and second excited 
+ states in Nd nuclei and, in the right panel, the isomer 



3 





40 


£ 

M— 




CM 

O 


36 




oz 


© 




CM 

1 


28 


"cn* 

+ 


24 






© 


20 



i — i — r 

o— o 



.6 



i — r 
(a) I 



o 



Nd 



J I I L 



\ ■ 

\ 



84 86 88 90 92 94 
N 



t — r 

o 
/ 



- / 

/ 

o 



t — i — i — r 

(b) 



o 



o. 



J I I I I L 



40 



- 30 



- 20 



- 10 



84 86 88 90 92 94 96 
N 



FIG. 2: (Color online) Calculated differences between squares 
of ground-state charge radii: {r%) Q +{A + 2) — (r^) +(A) (left 

panel), and isomer shifts (r;?) 2 + — (r%} + (right panel), as 
functions of neutron number in Nd isotopes. 
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FIG. 3: (Color online) Evolution of the first and second ex- 
cited + states (left panel), and the isomer shifts (rj?) Q + — 

(r c ) + (right panel) with neutron number in Nd isotopes. 
Microscopic values calculated with the PC-F1 energy density 
functional are compared to data [22l |. 
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(?" 2 ) +. The excitation energies of both 



and 0;j~ exhibit a pronounced dip at TV = 90, which can 
be attributed to the softness of the potential with respect 
to f3 deformation in 150 Nd. For lighter nuclei, i.e. toward 
spherical shapes, 0^ and 0^ display the structure of two- 
and three-phonon states, respectively. The axially de- 
formed 152 ' 154 > 156 Nd are characterized by strong prolate 
minima, stiffcr potentials, and the positions of (3 and (3(3 
bands arc shifted to higher energies. This behavior of the 
0^" state in Nd-Sm-Gd isotopes was also predicted in the 
phcnomcnological analysis of the transition between the 
vibrational, SU(5), and the rotational SU(3), limits of 
the interacting boson model (IBM) [2l|. We note that, 
with the exception of the very low 0^ state in 146 Nd, the 
calculated excitation energies E Q + are also in quantita- 
tive agreement with experimental values [22[ ■ 

The microscopic calculation predicts a very interest- 
ing evolution of the isomer shift (j" 2 ) + — {r^) +- After a 
steep rise with neutron number for TV < 90, the isomer 
shift actually changes sign between TV = 90 and TV = 92, 
i.e. in 152 . 154 < 156 Nd the charge radius of the 0^" state is 
smaller than that of the ground state. As function of 
the control parameter (number of neutrons), the isomer 
shift displays a behavior characteristic for a first-order 
phase transition [121 ]. Evidence for deformation crossing 
near the first-order shape-phase transition in 152 - 156 Gd 
has recently been reported in a study that used exper- 
imental B(E2) values to extract the model-independent 
quadrupole shape invariants, which provide a measure of 
the /3-dcformation 



finitencss of the nuclear system. In general, the charac- 
teristic behavior of order parameters at the point of QPT 
is more pronounced than in the case of Ising-type Hamil- 
tonians representing systems of interacting bosons, espe- 
cially for a realistic number of bosons, i.e. TVb ~ 5 — 10 
for medium heavy nuclei. In the latter case the dis- 
continuities arc smoothed out so that, qualitatively, a 
first-order phase transition might actually appear like a 
second-order one [l2[ ■ 

Shape transitions and change in radii are also reflected 
in the transition matrix elements of the electric monopolc 
T(E0) operator j24[. In the study of sharply rising E0 
strength in transitional nuclei [25| . based on a general 
interacting-boson model Hamiltonian of Ising-type, it has 
been shown that 0^ — > 0^ transitions provide a clear 
signature of phase transitional behavior in finite nuclei. 
In shape transition regions p 2 (75'0;0j — » 0*) displays a 
steep rise, and then remains large for well deformed nu- 
clei. Relative and absolute E0 transition strengths on 
the transitional path between the X(5) solution and the 
rigid rotor limit have recently been evaluated using the 
confined /3-soft (axially symmetric) rotor model (CBS) 
[26| . and it has been shown that absolute E0 transition 
strengths are reduced with increasing potential stiffness 
toward zero in the rigid rotor limit. 

The E0 operator can be expressed in terms of single- 

nucleon degrees of freedom as: T(E0) = J2 e k r ki where 

fc 

6^ is the charge of the fcth nucleon, and is its position 
relative to the center of mass of the nucleus. For the 



1|. Note that, although the results transition 2 -> 1 the absolute E0 strength is defined 



correspond to a realistic calculation of ground-states and 
collective excitation spectra of Nd nuclei, both isomer 
shifts, i.e. (r 2 ) 2 + - (r 2 ) +, and (r 2 ) + - (r 2 ) +, exhibit 
very sharp discontinuities at TV = 90. The first-order 
phase transition appears not to be smoothed out by the 



p 2 (E0;0+^Q+ 



(0+|T(E0)|0+ 



2i? 2 



(1) 



where R is the nuclear radius, R ~ 1 .2A 1 / 3 fm. Fig. 0] 
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FIG. 4: (Color online) The calculated monopole transition 
strength p 2 {E0; Of -> 0^ ) as a function of neutron number N 
in Nd isotopes. 



shows the calculated values p 2 (E0; Of — > 0^) as a func- 
tion of neutron number iV. Bare charges have been used 
in the calculation, i.e. e p = e and e„ = 0. The monopole 
transition strengths exhibits a markedly sharp increase 
toward the point of phase transition at N = 90, and 
the p 2 (E0;0+ -> 0+) values remain rather large in the 
well-deformed nuclei 152 454,156]\t ( j^ a resu i^ similar to that 
obtained in the schematic interacting-boson model calcu- 
lation of Ref. [25]. This behavior of p 2 (E0; 0+ 



characteristic for an order parameter at the point of first- 
order QPT. In terms of absolute values, one notes that 
even without introducing effective charges, the calculated 
p 2 (i?0;0^ — » Of) are in qualitative agreement with the 



available experimental values for Sm and Gd nuclei 24 1 . 

In conclusion, a microscopic calculation of observables 
related to order parameters for a first-order nuclear QPT 
between spherical and axially deformed shapes have been 
performed. Starting from self-consistent triaxial mean- 
field binding energy maps in the (3 — 7 plane for a se- 
quence of even-even Nd isotopes with neutron number 
TV = 84 — 96, a set of observables has been computed 
using collective wave functions obtained by diagonaliza- 
tion of the corresponding five-dimensional Hamiltonian 
for quadrupole vibrational and rotational degrees of free- 
dom. The energy gap between the ground state and the 
excited vibrational states with zero angular momentum, 
the isomer shifts {r 2 c ) 2 + - {r 2 ) Q + and {r 2 ) a + - (r 2 } +, 

and the monopole transition strengths p 2 {E0; Of — * 0+) 
exhibit pronounced discontinuities at N = 90, character- 
istic of first-order QPT. Even though the calculation has 
been carried out for a finite number of nucleons, the phase 
transition does not appear to be significantly smoothed 
out by the finiteness of the nuclear system. Together 
with the results reported in Refs. tne present anal- 

ysis has shown that the microscopic framework based 
on universal energy density functionals provides a fully 
consistent description of nuclear shape QPT in the rare- 
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